Systematic Field Theory of the RNA Glass Transition 
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We prove that the Lassig-Wiese (LW) field theory for the freezing transition of the secondary structure 
of random RNA is renormalizable to all orders in perturbation theory. The proof relies on a formulation 
of the model in terms of random walks and on the use of the multilocal operator product expansion. 
Renormalizability allows us to work in the simpler scheme of open polymers, and to obtain the critical 
exponents at 2-loop order. It also allows to prove some exact exponent identities, conjectured in LW. 
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Together with DNA and proteins, RNA plays a key role 
in biology. As such it is important to understand their spa- 
tial conformations. While for proteins the lowest-energy 
fold depends strongly on the chemical constitution, and is 
only tractable numerically, the problem for RNA is sim- 
pler, due to a clear separation in energy scale between pri- 
mary structure (the sequence), secondary structure (pairing 
of bases in a fold) and tertiary structure (embedding of a 
fold in 3-d space). The homopolymer problem (all bases 
identical) was solved in 1968 by de Gennes |1]. He finds 
that the pairing-probability of two RNA-bases with labels 
s and t, counted along the backbone, scales like V(s, t) ~ 
\s — t\~ Po , with Pq = 3/2. Real RNA molecules however 
consist of a sequence of 4 different bases, and their optimal 
fold depends on this sequence. Experimentally important 
(see e.g. ||2|]) is further the observation, that pairings (s,t) 
and (s',f) are either nested (s < s' < t' < t) or inde- 
pendent (s < t < s' < t'), which graphically amounts to 
the rule to draw the sequence and the pairings on the plane 
without self-intersections (planarity). While the problem 
of a biological RNA-sequence is best solved numerically, 
for reference it is crucial to understand the physics of (pla- 
nar) pairings of a random sequence. This was pioniered 
by Bundschuh and Hwa (BH) 1 3]. They consider a random 
pairing model with partition function (/3 = 1/kgT) 
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2^ rj(s,t)<P(s,t) 

\<s<t<L 
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which is defined as a sum over all planar pairings <£, such 
that <3?(s, t) = 1 if (s, t) is a Watson-Crick pair, and oth- 
erwise. The pair energy r](s, t) is considered as a quenched 
Gaussian disorder variable r](s, t), with 



rj(s,t) = f, 7](s,t)rj(u,v)—f 2 = a 2 5(s—u)5(t—v). 

Note that this is an additional approximation from the 
model of a random sequence [4]. A key feature of the 
above model is that there should be a continuous freez- 
ing transition between a weak-disorder phase, at large 



scales undistinguishable from the homopolymer case, and 
a strong-disorder glass phase with non-trivial scaling, and 
of possible biological relevance since the conformation and 
properties of RNA depends on the sequence disorder, i.e. 
on the primary structure. This glass phase appears in the 
BH solution of Q for the n = 2 replica case (instead of 
n = relevant for the disordered case) and in numerical 
studies at strong disorder Although the initial BH 

model is quite simple, this strong disorder phase of random 
RNA appears to be highly non-trivial and quite difficult to 
study, making it a challenging problem. 

In [6] Lassig and Wiese (LW) pioneered a field theoreti- 
cal (FT) approach for the transition to the glass phase. They 
showed their model to be renormalizable at first order in 
perturbation theory and calculated the critical exponents. 
Using a locking argument (see below), the scaling expo- 
nents for random RNA in the strong disorder phase were 
derived, in good agreement with numerics Si. 

It is important to understand if this approach defines a 
consistent theory to all orders, and if the estimates of @l 
for the scaling exponents are reliable. We achieve this goal 
here. Using a formulation of the LW model in terms of in- 
teracting random walks (RW) in d = 3 dimensions, and 
FT tools developped for polymers and membranes 0, 0] 
we show that the LW model is renormalizable to all orders. 
Our formulation is more convenient for calculations and al- 
lows us to derive new scaling relations between exponents, 
and to calculate critical exponents at second order. 

The FT model of LW is defined through perturbation the- 
ory in the disorder strength g = a 2 , the replica trick, and 
the continuum limit where L — > oo. One introduces n 
replicas of the system, labelled by a = 1, . . . , n. For the 
free model (no disorder, i.e. a = 0) the replicas are not 
coupled and the expectation value of a product of N con- 
tact operators & a (si,ti) can be computed exactly. It de- 
scribes the constrained configuration with iV fixed pairings 
(si, ti) (i = 1, • • • , N), i.e. N + 1 subrings of backbone 
length £ , ■ ■ ■ ,£ N (with L = £ + • • • + £ N ). As discussed 
above, this e.v. is 
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(* ( Sl> h) ■ ■ ■ $ a (s N , t N )) = £- pa l- po ■ ■ ■ l- po (2) 
with p Q = 3/2 if the (s,t)'s form a planar pairing, and 



2 



otherwise. The partition function for a single free RNA 



strand then is > = (1)„ = L~ Po . 

The average over the disorder 77 generates an attractive 
interaction between two replicas, 



-9o 

a<0 



Ks<t<L 



with coupling g oc a 2 and the overlap operator 

*a/s(s,t) = * Q (S|*)*/j(s,t) • 

The quenched disorder average is obtained for n 
averaged free energy F for a single strand is 

d 
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Z n is the partition function for n interacting replicas. Simi- 
larly, the average of an observable A, A, is the n — > limit 
of the average in the interacting theory. These observables 
are calculated as perturbative series in the disorder strength 
<7o- They suffer from short-distance UV divergences. Tak- 
ing po as an analytic regularization parameter, p < 3/2, 
LW show at first order that these UV divergences are poles 
in e = 2p — 2 at e = and that the theory is 1-loop renor- 
malizable at e = (p = 1). An UV- finite renormalized 
theory is defined through a renormalization of the disorder 
strength g and of the backbone length L. This allows to 
compute at first order the RG (3 function for the disorder 
strength g. It is found to have a UV fixed point g* > for 
e > 0, in particular for the physical case e = 1, n = 0. g* 
corresponds to the RNA freezing transition. LW compute 
also the scaling dimensions p* and 8* of the operators 
and \& at the transition. 

Our goal is to construct a FT which reproduces 0. For 
this we note that i~ Pa is the return probability at proper 
time £ for a free random walk (RW) in M d with d = 2p . 
Thus we introduce n independent RW's, r a . a = 1, . . . , n 
(r a (s) = {r£(s); p = 1, . . . , d}). In order to keep only 
planar pairings we use nxN pairs of auxiliary fields 7" (s) 
and 7"(s) (a = 1, . . . , N). The free model is given by the 
action ( ' := d/ds) 



a=l 
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n N r .L 



a — 1 a— 1 



(6) 



The propagators for r and 7, 7 are 



(7) 



<7 a Q (s)7f (*))o = ^^(* " «) , (77) = <77> = , 

where 8(s) = 1 if s > 0, and otherwise. 

The key point is that in the large- N limit, the observ- 
ables for a strand of length L in the LW model correspond 
to the partition function for a closed RW in our model, with 
specific boundary conditions (the end points are fixed and 



there is a creation operator 7 at the origin and an annihila- 
tion operator 7 at the end). The contact operator $ changes 
to 



N 

a,b 

The pair-contact operator fy a p(s,t) is still given by (0), 
and the interaction by The auxiliary fields 7 and 7 al- 
low to select planar diagrams by taking — ► 00. For the 
analysis of the UV-divergences, they are mere spectators. 
Their importance is that they allow to write an action, and 
thus to apply established tools to prove renormalizability, 
and to obtain exponents at higher orders. For the sake of 
simplicity, we shall not write the 7's explicitly in the fol- 
lowing. We also note that a 1/A-expansion is feasible, 
similar in spirit to |9|| for the homopolymer problem. 

The model defined by (0 belongs to a class of theo- 
ries with multilocal interactions, including the Edwards 
model of polymers and Self Avoiding Manifolds (SAM) 
0, 0| . Its short-distance singularities can be studied by 
the same Multilocal Operator Product Expansion (MOPE). 
Indeed, the operator ^ is a product of bi-local operators 
S[r a (s) — r a (t)] and of auxiliary fields 77. These auxil- 
iary fields have a very simple propagator and a trivial short- 
distance expansion, which is a product of 8 functions, mul- 
tiplying the MOPE of the S's. Let us give as examples the 
configurations which encode the UV singularities relevant 
at one loop. The short-distance behavior of a single ^ is 
given by 

■V al3 (u,v) ~ \u-v\~ d l (8) 
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and is depicted graphically as 



Similarly, two 'f's can coalesce into a single in three 
ways. Firstly as 

^ a p{u,v)^ a0 (u',v') ~ C(u',u;v',v)^f a/3 (u,v) H 

u'— *u,v f —*v 

(9) 

that we depict as 



with C(u',u;v',v) = (\u' — u\ + \v' — v \) the corre- 
sponding MOPE coefficient. Secondly as 

*ofl(«,tt')^(«y) - D(u,u',v)y af3 (v,v') + ■■■ (10) 



u — >v,u' — >v 



with D(u,u',v) 



u —u 



d if u<u'<v or v<u<u', and 



otherwise, that we depict as 



+ 
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Thirdly as 

■® M (u,u')^ a p(v,v') ~ E(u,u',v)^ af3 (v,v') + ■■■ (11) 

u — >v,u' — >v 

with E(u,u',v) =D(u,u',v), that we depict as 



+ 



The perturbative expansion involves expectation values of 
integrals of products of operators. The short-distance 
contribution u— m' for a single ^f(u, u') in these integrals 
is given by (JHJ and produces an UV divergence. The first 
term in l|8j gives a UV pole at d = 1 proportional to the 
insertion of the unit operator 1, while the second one gives 
a pole did = 2 proportional to the operator r 2 . Similarly, 
considering now the integrals involving two ^ operators, 
the integrals over u and v! (v and v' fixed) in (I9I10I1U 
give UV poles at d = 2, proportional to the operator 
In both cases, the subdominant terms (represented by the 
• • • ) in the MOPE involve higher dimensional multilocal 
operators, but do not give any UV pole at d < 2. Note 
that although the l.h.s of (lilt involves three replicas, the 
dominant term on the r.h.s. involves only the overlap op- 
erator between two replicas. Also note that (llOt does not 
contribute for polymers or SAM's, since there is a third 
non-planar diagram, and the sum of all cancel. 

This analysis of the UV divergences through the MOPE 
at first order gives the same results as in LW. It shows that 
our model is renormalizable to one loop at d = 2, as ex- 
pected from the existence of an action, and dimensional 
analysis. 

Our formulation and the MOPE allow to extend this 
analysis to all orders of perturbation theory fioll . The di- 
mension d = 2p of imbedding space is a dimensional reg- 
ularization parameter, and short-distance UV divergences 
appear always as poles in e = d — 2 = 2p — 2. We have 
shown that the theory is UV finite for e < (apart from 
a trivial "vacuum energy divergence" proportional to the 
unity operator 1). For e = the only UV divergences are 
proportional to the local operator r 2 and to the bilocal op- 
erator ^> a j3 (a 7^ (3). The MOPE also generates multilocal 
operators involving more than 2 replicas, for instance the 
3-replica operator II Q/ 3 7 = $ ct $ / 3$ 7 . However these oper- 
ators are not associated to UV divergences, and correspond 
to irrelevant couplings. The crucial point in this analysis is 
that, since the interaction in the model involves 2 different 
replicas a ^ (3, no UV divergence appears which is pro- 
portional to the single -replica operator 3> a (u, v), although 
$ is a "dangerously" marginal operator (it is marginal at 
e = and relevant as soon as e > 0, like 

The renormalized UV finite theory is defined through the 
renormalized action S R 



E 



J0 4 a<P Jo<u<v<L 



Z and Z g are the wave -function and coupling-constant 
counterterms, and are series in g R whose coefficients con- 
tain the poles in 1/e. p is the renormalization mass scale. 
The renormalized field r and coupling g R are related to 
the bare ones r B and g B by r B = Z x I 2 y and g B = 
g R Z g Z d p~ e (this differs from the LW scheme where s is 
renormalized instead of r and where g' B = g R Z g Z 2 p~ £ ). 
The RG (3 function for the coupling (3 g and the scaling di- 
mension Xr f° r tne field r in length-units are 

dg R 



0y 



dp 



1 / dlogZ 
^=2( 1 + ^^t 



(13) 



Since we proved renormalizability, and identified all pos- 
sible counterterms, we can simplify calculations by us- 
ing open RWs instead of closed ones, eliminating the 5- 
function for the closure. Although not all correlation func- 
tions are directly interpre table in terms of RNA strands, 
they are renormalized by the same counterterms, except 
for one additional boundary term % x f° r eac h end of the 
RW, S° R pen = S c R loscd + 2Zl The (Fourier transformed) 
partition function for a single free open RW then is 



= (lie 



-nq L 



a=l 



(14) 

The first-order correction to (1141 in the disorder strength is 
given by the following diagram 

n(n — 1 
9o ^ — 



Z 



e <j (2\v-u\-nL)\_\-d ^ 



0<u<v<L 



(15) 

The last integral is UV divergent when u — > v (bulk) and 
u, v — > or u, v — > L (boundary). The corresponding UV 
pole in e = has residue (1 — 2Lq 2 ) exp(—nq 2 L). This 
partition function for a single open RW is renormalized as 

Z%\q,g R ) = Z^iZ-^g^e-^ . (16) 

This implies that in the MS scheme the counterterms are at 
first order Z = 1 + g R (n — l)/e, = g R (n — l)/4e. 

To compute Z g it is simpler to consider the (connected) 
partition function Z^ for two distinct interacting open 
RWs. At first order in g it is given by 
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9o 



(17) 



ZL 



At order g 2 there are 4 UV divergent diagrams, with MOPE 
given in JTOl . and (fill . 



(12) 



and the corresponding function is renormalized as 

Z R 2 \g R ) = Z- d Z ( i\g B )e- 4 ^ . (18) 

Care was taken to account for the (missing) zero-mode due 
to the ^-distribution between the 2 replicas, resulting in the 
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FIG. 1: Results for 8* (black) and ip* (grey) at 1- (dashed) and 
2-loop (solid) order. 



factor of 7Lr d . To subtract the UV pole at e = the coun- 
terterm is Z ff = 1 + gu(7 — 4n)/e. 

This scheme can be continued to 2-loop order lHoll . To 
simplify the results we subtract minimally Z and Z s Z d . 
Also p(g r ), the dimension of <£, can be calculated by con- 
sidering a 2-RW partition function Z$ with one <I> connect- 
ing the 2 RW's. We obtain the RG functions at 2 loops: 
-eg R + (5-2n)g 2 R + (3-2n)(5 

1 n — 1 (n 

2 " ~^ 9R " 



Xr(gn) 



l)(4-3n) 



2n)^ 
(19) 



1 + - + {n 



(n-l)(3 



At one loop our results agree with those of LW |6], upon 
identifying /3 = /3 LW , and Xr = 1/(27 LW ). In the physical 
case of random RNA (n = 0), our 2-loop results confirm 
the existence of a UV fixed point (FP) (in our scheme at 
g* = |e — ^e 2 ), describing the freezing transition. The 
anomalous dimensions of ^ and <I> at this FP are 

6* = A4g*) = 2 + P' g (g*) = 2-e-^e 2 (20) 

p * = A*(g*) = l + ^e + -^e 2 . (21) 

On Fig. [2 we plot 6?* (black) and p* (grey) from e = 
to e = 1 (physical case). The full curves are our 2-loop 
results, the dashed ones the 1-loop estimates of [6]. The 
2-loop corrections do not change much the estimate for p* , 
but do change it quantitatively for 8*. This defines another 
dependent exponent, the roughness £ of the height-field 
Kr) ■= J2 s< rJ2 t>r Hs,t) ~ rC, as C = 2 - p fl. 
Since $ a (s,t) < 1, (& a p) (s,t) < ($„(s,t)) and thus 
p < 6. This bound is clearly violated by our results for 
e > e c ~ 0.59. According to LW in this regime the two 
replicas are locked into a single conformation, p* = 6* and 
the exponents in the glass phase equal those at the transi- 
tion 

Pglass = P*|e=l ^ 1-36, Cglass = 2- p* ~ 0.64 . (22) 

Finding a small e 2 correction to p* is important, since it 
validates the estimates of |6] for the exponents of random 
RNA. Numerical results obtained by Krzakala et al. j5j] in 
agreement with Bundschuh et al. loL fllll give 

Pgiass * 1-34 ± 0.003, Cgiass 0.67 ± 0.02 (23) 



In it was also conjectured that the dimensions of and 
r are not independent. Our formalism shows that this is 
correct and gives an exact relation between \r an d P- We 
remark that the partition function Z$ for one $ connecting 
two strands is equivalent to that of two single strands, upon 
marking a single point on each strand, i.e. 

1 



Z^, — — 



q=0. 



(24) 



(25) 



At g*, together with Q + p = 2 [6], this gives 

C* = 2 - p* = (2 - e) X * r ■ 

In conclusion, our results for the RNA freezing transition 
are as follows: First we give a new field theoretical for- 
mulation of Lassig-Wiese [6], and prove that their model 
is renormalizable to all orders in perturbation theory. As a 
consequence we show that the e-expansion scheme of 0| 
is well defined. Second our formulation allows to simplify 
the perturbative calculations, in particular by considering 
open interacting RW's instead of closed RNA strands. We 
perform the first 2-loop calculation for the critical expo- 
nents 6 and p, and show that it does not much correct the 
LW estimate for p. Third, we derive a new scaling law re- 
lating the dimensions of the height field h and r. Finally 
let us mention that we have applied our methods to the de- 
naturation of random RNA under tension, which allows to 
calculate the extension force curve IH2ll . 

K.W. thanks M. Lassig for the stimulating collaboration 
which raised many of the questions addressed above. 
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